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Abstract 



■ Current solar and atmospheric neutrino oscillation data seem to favor a bimax- 

Ph. imal pattern for neutrino mixings where the matrix elements U e 2 and are of 

Q_i! order one, while U e s is much smaller. We show that such a pattern can be obtained 

quite easily in theories with "lopsided" mass matrices for the charged leptons and 
the down type quarks. A relation connecting the solar and atmospheric neutrino 
mixing angles is derived, tan 2 Q a tm — 1 + tan 2 so i, which predicts sin 2 26 atm ~ 0.97 
^ ■ corresponding to the best fit LMA solution for solar neutrinos. Predictive schemes 

in 50(10) realizing these ideas are presented. A new class of S'O(IO) models with 
lopsided mass matrices is found which makes use of an adjoint VEV along the I3R 
direction, rather than the traditional B — L direction. 



§1. Introduction 



Recent data seems somewhat to favor either the LMA (large mixing angle MSW) or 
the LOW solution to the solar neutrino problem over the small mixing angle solution 
P], Taken together with the atmospheric neutrino results B and the CHOOZ reactor 
experiment Q, this would imply a so-called "bimaximal" pattern of mixing, with U e 2 and 

large and U e s small (U is the leptonic mixing matrix) ||. Another possibly significant 
feature of the data is that the atmospheric neutrino mixing angle is not merely large, 
but seems to be nearly maximal. The best-fit value at present is sin 2 29 atm - 1.0 § 

In this paper we make several points relevant to these observations. 

1. The so-called "lopsided" models f| [7|, §, ||, [L0], [TTJ] provide a very simple way of 



accounting for the bimaximal pattern of mixing, and in particular have no diffi- 
culty in obtaining the LMA solar solution, unlike certain other kinds of bimaximal 
schemes. 

By combining the idea of lopsided mass matrices with a nonabelian flavor symme- 
try one can explain in a simple way the near maximality of 6 atm . In the simplest 
case one obtains the relation tan 2 8 atm = 1 + tan 2 (9 so ;. In the limit of small solar 



angle this gives maximal atmospheric angle; while for the best-fit LMA value |T2 
of tan 2 6 so i = 0.4, it gives sin 2 26 atm = 0.97. 

The lopsided bimaximal idea can be straightforwardly implemented in the context 
of 50(10), and in that case a prediction relating quark masses and mixings to the 
atmospheric neutrino mixing angle arises. 

A new class of predictive S*O(10) models for quark and lepton masses is found 
which makes use of an adjoint VEV along the I^r direction. Bimaximal mixing 
pattern for neutrinos can be obtained easily in this class of models, along with 
several predictions relating the charged fermion masses and mixings. This provides 
a new way of looking at quark and lepton masses in S'O(IO), different from the 
traditional way where the VEV of the adjoint points along the B — L direction. 



§2. Bimaximal mixing 

Imagine that the leptonic mixing angles come primarily from the diagonalization of 
the charged lepton mass matrix L, which has the following "lopsided" form: 

L= e \m D . (1) 

\P' P 1 J 

Here p' ~ p ~ 1 , whereas e < 1. The dashes represent elements that are small compared 
to the ones shown. The convention being used is that the left-handed lepton fields 
multiply the mass matrix from the right. The diagonalization of this matrix can be 
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done in stages, the first stage being to rotate in the space of £ 2 and l x by an angle 
which we will call 9 S , satisfying tan9 s = p'/p. This brings the matrix to the form 



u 



a 



(2) 



where a = \J p' 2 + p' 2 . (Note that all parameters shown in L can be made real by field 
redefinitions.) The next stage is to rotate in the space of £3 and the new £ 2 °y an angle 
which we will call 9 a , satisfying tan# a = a. This brings the matrix to the form 



/ - 



V 



1 







m D . 



(3) 



The rotations needed to complete the diagonalization involve only small rotations 
of the left-handed leptons, and we will therefore neglect them. (An important point is 
that the (2,1) and (2,2) elements of L were assumed small compared to e. Otherwise, 
there would still be required a large rotation in the 1-2 plane to diagonalize L", and that 
would induce a large U e ^.) The unitary matrix Ul required to diagonalize L^L is thus 
approximately 



10 

cos 9 a - sin 9 Q 
sin 9 a cos 9 a 















s 






C a 



(4) 



The full leptonic mixing matrix is given by Umns = U L U U , where U v is the unitary 
matrix required to diagonalize the neutrino mass matrix. However, since we are assuming 
U v = I, Umns is given approximately by the matrix in Eq. (4), which has the bimaximal 
mixing pattern, with U e2 and U^ both of order unity and U e % small. A very important 
point is that no constraint whatsoever has had to be placed on the neutrino masses. 
The questions of neutrino mixing and neutrino mass are completely decoupled in this 
scenario. This means, in particular, that there is no difficulty in obtaining the neutrino 
mass ratios appropriate to any of the large angle solar solutions, LMA, LOW, and VAC. 

In many published models the bimaximal mixing comes from the diagonalization of 
the neutrino mass matrix fll3| . It is instructive to compare the present idea to some of 
these other approaches. Consider the following three forms of M u , the light neutrino 
mass matrix obtained after seesaw diagonalization. 



M A ~ 



a 
a 



a 
1 



P' 2 PP' P' 
PP' P 2 P 
P' P 1 



M° ~ 

±VJ -y 



- 1 

1 - 
o — 



a 



(5) 
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As before, the dashes indicate elements smaller than the ones explicitly shown, and p, 
p' and a are assumed to be of order unity. 

In matrix M^, a large rotation angle, satisfying tan# a = a, is required to diagonalize 
the 2-3 block. This produces a large atmospheric neutrino mixing angle. The magnitude 
of the solar neutrino mixing angle depends on the magnitude of the small elements in 
M^, and may also be large. Because of the approximately "factorized" or rank-1 form of 
the 2-3 block of this matrix, there is only one large mass eigenvalue, so that the desired 
hierarchy 777.1,777.2 <C results. Matrix can thus give a satisfactory bimaximal 
mixing. However, there is a price to be paid for this: in order for both the mixing 
angles and the neutrino masses to come out right a certain precise relationship had to 
be assumed to exist among the elements of M v — namely, the approximately factorized 
or rank-1 structure of the 2-3 block. Moreover, for the solar neutrino mixing angle to 
be large, further assumptions have to be made about the small elements of M„. 

Matrix M„ has the apparent advantage over matrix that both the atmospheric 
and the solar neutrino mixing angles automatically come out to be order one. However, 
it does not give a realistic bimaximal scheme, since U e s is of order one rather than 
small. The reason is the following. To diagonalize Mjf , requires first rotating in the 
1-2 plane by an angle with tan0 s = p'/p, and then in the 2-3 plane by an angle with 
tan 9 a = o = Vp' 2 + P 1 ■ This is the same as what was required to diagonalize the charged 
lepton mass matrix in Eq. (1). However, because one is diagonalizing the neutrino mass 
matrix in this case, the resulting MNS matrix is the adjoint of what was obtained in 
Eq. (4). (Recall that Umns — U\jJ u .) Thus, here U e ^ = sin 9 a sin 9 S . Moreover, as in 
the previous example, matrix M„ requires a form in which the elements are in a special 
precise relationship to each other. 

Matrix is the typical "inverted hierarchy" form |14]], with an approximate L e — 

— L T symmetry and automatically gives bimaximal mixing, with C/ e 3 very small. This 
can be seen as follows. The matrix can be diagonalized in stages, as in the other 
examples. In this case, however, the first stage is to rotate in the 2-3 plane by an angle 
6 a such that tan 9 a = a. This brings the matrix to a "pseudo-Dirac" form with large and 
equal (1,2) and (2,1) elements and all other elements small. The next stage is a rotation 
by 7r/4 in the 1-2 plane. Thus, U u has the form given in Eq. (4) with 9 S = tt/4. The 
fact that the solar neutrino angle typically comes out very close to maximal is certainly 
acceptable for the LOW and VAC solutions, but may not be acceptable for the LMA 
solution if it turns out that the LMA fits require tan 2 9 so i to be significantly smaller than 
one. At the moment, the best-fit LMA value is tan 2 9 so i ~ 0.4, but maximal mixing is 



within the 99% confidence level contours given in [12 



One sees from these comparisons that obtaining bimaximal mixing from the diago- 
nalization of the charged lepton mass matrix simplifies the problem by neatly separating 
the questions of neutrino mass and neutrino mixing. 

§3. Nearly maximal atmospheric neutrino mixing 

At present the best fit to the atmospheric neutrino angle is sin 2 29 a t m ~ 1.0. This is 



3 



difficult to obtain as a prediction from models. Several types of model, indeed, predict 
that the solar angle should be very close to maximal — for instance, the "inverted 
hierarchy" models [FJ[] just described and "flavor democracy" models [D|. But maximal 



atmospheric neutrino mixing is much harder to achieve. The reason is simple. The 
most obvious way to get nearly maximal mixing of two neutrino flavors is by a pseudo- 

Dirac form of the mass matrix: ( ^ ^ V with i5,5' < 1. This form also gives nearly 

degenerate neutrinos. Therefore, if such a form is assumed for the 1-2 block of M v , 
to give maximal solar neutrino mixing angle, it also typically gives Am^ <C Am^ m , 
as desired. However, if the 2-3 block of the neutrino mass matrix is assumed to have 
a pseudo-Dirac form, to give maximal atmospheric neutrino mixing angle, it typically 
gives Aml tm <C Am^,, which is wrong. 

It is quite difficult to find a form of M v that both gives maximal atmospheric neutrino 
mixing angle and Am^ <C Am^ (m . The ingenious model of Ref. Jl6| shows what is 
required to obtain this result. 

It is much easier to obtain maximal atmospheric neutrino mixing through the charged 
lepton mass matrix [17|], precisely because that decouples the neutrino mixing pattern 
from the neutrino mass pattern. Consequently, the requirement that Am^ <C Am 2 tm 
presents no difficulty. All that is needed is that there be a nonabelian symmetry relating 
pT and r~ so that the parameter p in Eq. (1) comes out to have magnitude 1. This can 
be done in various ways. One possibility is that (px,r£) = ipl form a doublet of the 
permutation group 5*3, while the e~l is a singlet. If the S3 is broken by a doublet "flavon" 
field Xi with its VEV given by (xi) = (1, i) (this form of the VEV can emerge from certain 
simple forms of the Higgs potential as shown below), then the desired (3,2) and (3,3) 
elements of the matrix L given in Eq. (1) can arise from the term Tj^i^J (xi) {Hd) . This 
will make \p\ = 1 in Eq. (1). 

Simple Higgs potentials can be constructed with flavon fields that have the desired 
VEV structure. As an example, consider the potential corresponding to an S 3 doublet 
flavons x = (Xi>X2)- The renormalizable potential for x that is invariant under S3 as 
well U(l) symmetry is [113] 



V(x) = f^ixlXi + X2X2) + \X1X1 + X*2X2) 2 + A 2 (xIX2 -xlxif 

+ H(xlX2 + xlxif + (xtxi - X2X2) 2 ] • (6) 

The VEVs can be parametrized as (xi) = rcos#, (xz) = rsin0e 1 ^. Minimization of V 
with respect to and 9 leads to <fi = ±tt/2 and 9 = ±7r/4, corresponding to (A 2 + A 3 ) 
having positive sign. This is the desired VEV, written as (x) = r(l, i). Realistic charged 
fermion masses can be induced by making use of flavon fileds which get VEVsof the form 
(1, 0) and (0, 1) in the space of the second and the third families. For an analysis of 
alternative ways of inducing this VEV structure in the context of supersymmetric models 



see Ref. 17 



If |p| = 1 in Eq. (1), the following relations obtain: tan 9 so i = p' / p = p', and 
tan 9 atm = a = V 1 p 1 + p' 2 = y/1 + p' 2 . Together these imply that 
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tan 2 9 atm = 1 + tan 2 6 soh (7) 
or, for the more usually quoted quantity, 

2 1 4" tan 2 9 so i 

sin 26 atm = , — y~ — - ? . (8) 



:i + |tan 2 9 sol/ 



One sees that as tan 2 6 so i varies between and 1, sin 2 26' a t m varies between 1 and 
8/9. (The point tan 2 9 so i = 1, sin 2 29 atm = 8/9, is the same as the prediction of flavor 
democracy models in the pure flavor democracy limit.) For the currently favored "best 
fit" value of tan 2 6 so i ~ 0.4, we have sin 2 29 atm ~ 0.97, in excellent agreement with data. 

§4. Embedding in grand unified models 

One of the main virtues of lopsided mass matrices, which has been emphasized in 
the literature || ||, |9|, PH , is that in the context of grand unified theories they very 
elegantly account for the disparity between the observed 2-3 mixings in the quark and 
lepton sectors, i.e. the fact that = 0.7 whereas V c b = 0.04. The explanation lies in 
the fact that SU(5) relates the charged lepton mass matrix L to the transpose of the 
down quark mass matrix D. In fact, in the "minimal" SU(5) model L = D T exactly. In 
the form of L shown in Eq. (1), it is the 0(1) elements p and p' that control the mixing 
of the left-handed fields and give large U^. However, if D is similar to the transpose 
of this form, then it is the small entry e (cf. Eq. (1)) that controls the mixing of the 
left-handed down quarks of the second and third family, namely V c b- It should be noted 
that SU(5) relates L only to D, and not to the up quark mass matrix U or the neutrino 
Dirac mass matrix N. Therefore, one expects that D should be lopsided if L is, but 
there is no reason to suppose that U and N are. In fact, in lopsided models that give 
a good account of quark and lepton masses and mixings, only D and L are assumed to 
have lopsided forms. This is true also of the realistic S'O(IO) lopsided models that have 
been constructed || [10 . 



Where the form in Eq. (1) differs from most published lopsided models is in the large 
element p 1 . (However, the model of Ref. 0, had a form for L much like Eq. (1), with 
an entire row of large elements.) The presence of the large element p' puts significant 
constraints on the building of realistic models of quark and lepton masses. The point 
has to do with the so-called Georgi-Jarlskog factors of 3: m s Pd and rrid ~ 3m e (at 

the GUT scale) |19[]. If the charged lepton mass matrix has the form shown in Eq. (1), 
then the simplest way to get the first Georgi-Jarlskog factor is by assuming that D has 
the form: 




D = - p \m D . (9) 

V - -e/3 1 J 

The factor of —1/3 in the e term relative to the corresponding term in L is easily 
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explained as being due to the 5*0(10) generator B — L. Exactly this factor appears in 
the 50(10) lopsided models of Refs. 



ID | 



To get the second Georgi-Jarlskog factor of three requires, as is well-known, that 
det D = det L. Barring some accidental cancellations, this forces the (1,1) and (1,2) 
elements of L (and correspondingly the (1,1) and (2,1) elements of D) to vanish, or at 
least to be negligibly small. One can also, without loss of generality rotate to make the 
(2,1) element of L and the (1,2) element of D vanish. This leads to the virtually unique 
forms for D and L. It is straightforward to generalize the 50(10) model of Refs. || 10 
to obtain the following realistic mass matrices: 



N 



( 


U 

t ° 






5 



5' 



m D , 








K 



mu, 





( 







D = 





5 


P + e/3 




u 


-e/3 


l + « / 




( 







H 





- 


-e/3 mu- 




vo 


e/3 


1 / 



(10) 



These mass matrices arise from the following Yukawa terms: The entries denoted '1' 
come from (163163)10^. The 0(1) elements k, p, and p' come from (163l6/f)(16il6#), 
where i = 1,2,3, the multiplets in the parentheses are contracted into 5O(10) vectors; 
16h breaks 5O(10) down to SU(5), and the 16' H breaks the electroweak interactions. 
The elements e, which appear antisymmetrically, come from 16216310^45^, where the 
VEV of the adjoint Higgs lies in the B — L direction. Such an adjoint VEV is what would 
be desired to achieve doublet-triplet splitting without fine-tuning via the Dimopoulos- 
Wilczek mechanism in 5O(10) [p0[| . The elements 5 and 5' come from terms of the same 
form as the k, p and p' terms, but with of course different family indices. These are 
exactly the same kinds of operators that appear in the models of Refs. 0, [8|, [10[ . 

These mass matrices give a quite satisfactory fit to all the quark and lepton masses 
and mixings. In the approximation l~p~p'~/t^e^5~5', so that the observed 
mass hierarchy is correctly reproduced, the following mass relations are obtained at the 
GUT scale: 



nib = m T , 



m M /3, m d = 3m e , m u /m t = 0. 



The first three are the Georgi-Jarlskog relations, all of which work quite well when 
compared with experimental values of the masses. m u /m t is predicted to be zero by 
these forms. Experimentally, it is about 10 -5 , which is about two orders of magnitude 
less than the corresponding ratio for the down quarks, m^/mfe. A tiny non-zero value of 
m u can easily arise from some higher-dimension operator. 

The parameter k in Eq. (10) is necessary in order to have adequate CP violation in 
the CKM matrix. We may redefine 1 + k to be simply 1 with an appropriate redefinition 
of mo in Eq. (10). The parameter e in L and D of Eq. (10) will then be different from 
e in N and U. Let us then rename e appearing in L and D as ze. In this redefined 
notation (we denote the redefined 5,8', p, p' by the same symbols) we have rris/m^ = 
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a/(l + a 2 )(ze/3), V us 5'/(ze/3), V ub 5'/(l + a 2 ), V cb (e/3)(z/(l + a 2 ) - 1). From 
these relations, we obtain the following prediction: 

- tan29c(m s /m b ) 
tan# aim — — — 

^| Vub\ 

where 9c is the Cabibbo angle, and Q a tm is the atmospheric neutrino mixing angle that 
comes from the charged lepton matrix. (The contribution from the neutrino sector 
is assumed to be small.) Note that 9 atm is of order unity, as needed for atmospheric 
neutrino oscillations. 

The mixing parameter U e3 is predicted to be 

\U e3 \^ S -^\U^\^ (0.04 - 0.05) (13) 

where the factor sin 9c/ 3 arises from the small rotation needed to complete the diago- 
nalization of L" of Eq. (3). (The factor 3 arises because = 3m s .) This prediction 
will provide a test of this class of models. 

If the parameter z were equal to 1 (which will be the case when the entry k is absent 
in Eq. (10)), then there will be not enough CKM type CP violation in this model, as all 
the mixing angles become approximately real. This is true even when we allow for the 
parameter 8 to be complex, since there is a cancellation between the up and the down 
quark contribution in the phase of the CKM matrix. Allowing for z ^ 1 (or k ^ 0) 
leads to the desired CP violation, since z is complex. It is interesting to note that if 
z were equal to 1, the charm mass will be predicted to be m c (m c ) = (1.1 — 1.2) GeV 
M. Furthermore, the relation \V ub \ = (rn s /m b ) 2 \V us /V cb \ will follow, which is in good 
agreement with experimental values. 

§5. New class of lopsided mass matrices from an Isr adjoint 

In the preceding example we made essential use of an 5O(10) adjoint VEV along the 
B — L direction. Now we show that quite simple and predictive mass matrices can be 
derived in 5O(10) if the VEV of the single adjoint present in the model points along the 
direction (I 3R stands for the third component of the right-handed isospin). There 
is a simple and elegant realization of the lopsidedness of D and L in this scheme. It is 
worth noting that a single adjoint scalar with its VEV along I 3 r direction can lead to a 
natural doublet-triplet splitting, just as in the case of a single B — L adjoint plf . The 
ha adjoint also has some advantages is suppressing Higgsino-mediated proton decay in 
supersymmetric 50(10) [p2| . 

Consider the case where the B—L adjoint that was involved in generating the lopsided 
mass matrices of Eq. (10) is replaced by an I 3 r adjoint. The mass matrices will then 
have the form 
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(12) 



N = 



( 


u 

t 







5 



5' \ 



( 



m D , 



D = 



1 J 










5 

—e 



m D , 



\ 



m v , 



U = 



( 









\ 



(14) 



/ 



As in the previous section, the '1' entries arise from I63I63IO// coupling. There are 
two crucial differences compared to Eq. (10). The entry resulting from 16216310^45^ 
has now two group contractions. These two are denoted in Eq. (14) as e and e'. These 
entries are proportional to the I^r charge, so that they are equal in D and L (similarly in 
iV and U). Lopsided nature arises from the p entry generated through 16216 3 16h16'# 
coupling. The parameters p, p' are assumed to be much larger than e, e'. This model 
then predicts the following relatins: 



m b = m T , 



m s ^m^ m d m s m b = m d m^m T 



(15) 



The inequality for m s follows since m s /m b 



|ecr/(l + a 2 )\, while m^/m T = \e'a/(l + 
a 2 )\, where a = \J p 1 + p n . Thus, although the entries e, e', proportional to I^r do not 
distinguish L from D, and the p-type entries also by themselves do not distinguish L and 
D (these entries do not break 5*^7(5)), a combination of the two leads to the breaking of 
= m s relation, as desired. Unlike in Eq. (10), there is sufficient CP violation in the 
CKM matrix in this model even without an entry like k of Eq. (10). 

Working in the approximation 1 ~ p ~ p' ^> e ~ e' ^> 5 ~ 5', we obtain the following 
relations for the masses: m b = m T = VT+ a 2 m D , m s /m b = [(ae + S*p/a)\]/(l + a 2 ), 
m^/rrir = [\(ae' + 5*p/a)\/(l + a 2 )}, m d m s m b = m e m^m T , m c /m t = ee', m u /m t = 0, all 
at the unification scale. The CKM mixing angles are given by \V US \ = S'/(e + 5p/cr 2 ), 
\V ub \ ^5'/(l + a 2 ), \V cb \ = |e(2 + cr 2 )-5*p|/(l + o- 2 ). Here all parameters have been made 
real by field redefinitions, except 5. The rephasing invariant CP violation parameter r\ is 
given by r] = lm{(V ub V cs ) / (V us V cb )} 2eplm(5)/{a 2 (l+a 2 )\V cb \ 2 }. From these relations, 
we obtain the following prediction for the atmospheric neutrino oscillations: 



tan# atm — a — 



{m s /m h )\V us \/\V uh \ 



(16) 



which is analogous to Eq. (12). We also have a quantitative prediction for tan# so ; = p' / p. 
This can be seen by noting that a is determined from Eq. (16), e' from m M /m r , e from 
m c /m t , plm(<5) from 77, pReS from m s /m b , and 5' from V ub . The determinant relation 
m e m M m T = \pp'S\ then fixes p' / p. 

Consider the input parameters taking the following values. \V US \ ~ 0.215, \V ub \ ~ 
0.0036, \V cb \ ~ 0.0037, r] ~ 0.33, and m c (m c ) ~ 1.35 GeV, m t = 175 GeV, all at the 
weak scale, along with m s /m b \GUT — 1/65 (corresponding to m s (\ GeV) = 130 MeV). 
This gives the prediction t&n9 atrn = a ~ 1.034, t&n9 sol = 0.43 and U e3 ~ 0.06, all of 
which are in reasonably good agreement with current neutrino oscillation data. 
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In summary, we have presented simple realizations of bimaximal neutrino mixing 
pattern, making use of lopsided mass matrices for the fermions. This idea has a natural 
embedding in unified SO (10) models. We have presented two different realizations within 
50(10), one making use of the traditional B — L adjoint VEV, and a new class of models 
making use of an I 3 r adjoint VEV. We were also able to derive in a simple way the near 
maximal mixing for atmospheric neutrino oscillation angle, as given in Eq. (8). 
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